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ABSTRACT : 

This report presents the aerodynamic calculations and describes the 
mechanical features of a test rig that has been designed at the Turbo- 
propulsion Laboratory, Department of Aeronautics, for research work on a 
special type of centrifugal compressor. This so-called Hybrid compressor 
consists of a centrifugal rotor with l80° flow deflection in the meridio- 
nal plane followed by an axial flow diffusor. The report establishes a 
method to predict the off-design performance of the subject compressor, 
which can be applied also to conventional centrifugal compressors. 
Extensive use is made of modern programmable calculators, and programs 
are presented to show the effectiveness of these tools in engineering 
endeavors . 
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INTRODUCTION 



The idea of arranging an axial diffusor after a special centrifugal 
compressor rotor that turns the flow by l80° instead of 90°, as in 
conventional machines, has been conceived by Mr, George Derderian, Naval 
Air Systems Command, Code AIR-53622B. For this so-called Hybrid com- 
pressor concept, he was awarded U. S. Patent 3? 365? 892 dated January 30, 
1968. 

To verify the conditions that exist after a Hybrid impeller, and to 
investigate whether it is possible to convert the kinetic energy of the 
flow after such a rotor into useful pressure rise with acceptable effic- 
iency, a research project was started at the Turbopropulsion Laboratory, 
NPS , which was funded by the Naval Air Systems Command, under: AIRTASK 

A3303300/i86b/3F4i432301 . 

The present report gives the particulars of the design that led to 
the manufacture of the test rig which is now operative and being used to 
carry out the afore-mentioned research work. 

The project was supervised by Mr. K. H. Guttmann, AIR-330C, who has 
been extremely helpful in all phases of the program. The author is 
deeply grateful to Mr. Guttmann for the patience, understanding and 
considerateness which he received during a lengthy period of personal 
misfortune which delayed the research program and the completion of this 
report. 
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I. AERODYNAMIC DESIGN 



1. Impeller 

Because of its peculiar design the rotor of the Hybrid compressor of 
Fig. 1 is highly stressed. (See Section V.l) Therefore, it must have 
meridional blades at the discharge to eliminate additional blade bending 
stresses. Ahead of the rotor the absolute inlet velocity will be axial 
everywhere. Hence it is possible to establish the rotor dimensions with 
the method of Ref. 1 which is enclosed as Appendix A. The corresponding 
dimensions of the Hybrid rotor are shown in Fig. 1. 

The symbols used in the following are identical with those listed in 
Appendix A. For easier reference, the figures and equations of Appendix 

A are denoted by Fig. A .... and Eq. A Thus, Fig. A (3) is figure 

3, and Eq. A II (4) is equation II (l4) of Table II of Appendix A. 

A parametric study indicated that optimum conditions were obtained 
for the following design values (see Fig. l): 

hx> = 55 ° 

a, = 65° 

R l 0 / R 2 = °- 55 

R li/ R 2 =0-18 

The relatively small angle 9 1q of 55° is necessary to be able to limit 

the ratio R, /R„ to 0.55 and to obtain a sufficiently high flow rate so 
lo d. 

that the diffusor blade height b 2 does not become too small. For values 
of R 1 q /R 2 larger than 0.55 the meridional flow channel would also have 
very large curvatures at the blade tip contour of the Hybrid rotor. This 
would produce large radial velocity gradients at the rotor discharge. 

The choice of = 65° was predicated by the permissible loading of the 
axial diffusor blade rows, and the large losses that occur if is 
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larger than 65°. A radius ratio = 0.l8 of the hub at the inducer inlet 

was chosen to be able to arrange 17 rotor blades with a thickness of 
about 0.125 inches. With smaller ratios R-^/R 2 ^ no ^ P oss ible to 
have 17 blades and it would then be necessary to design the impeller 
with splitter vanes, a feature which did not seem to be warranted for the 
present tests. 

Preliminary investigation also showed that the peripheral rotor 
speed at the mean discharge radius R^ of Fig. 1 should not exceed about 
800 ft/s because of critical speed considerations if the rotor is made of 
aluminum. 

For air with v = c /c = 1.4 and a gas constant R_ = 53-35 (ft -lb)/ 

' p v G 

(lbm, °R)> the velocity of sound a Q at an assumed total inlet temperature 
T Q = 520°R is 

a Q = y R q T q = 1117.85 ft/s 



giving a ratio 

U 2 /a Q = 0.7156 



From Eq. A Il(4) the Mach number M^ of the relative velocity at the 

outer radius R n of the rotor inlet eye is then 
lo 






. 2 
sin 6 



lo 



^10 v 
R 2 a 0‘ 



Y - 1 2 

~2— cos ^ 



= 0.4842 



(l) 



For a weight flow rate w = m g (lbm/s) there is from Eq. A Il(9) 

2 - - (11.7944) (2) 



ttR. 



*B1 



0 



Figure A(3a) shows that the so-called slip factor p, depends of the 
number of rotor blades the ratio R^ q /R 2 r the rotor efficiency T^, 
and the flow coefficient cp-,. Minor corrections for p, must be applied if 
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the tip clearance ratio s/b^ differs from 0.05 (see Fig. A (3b)) and if 
the flow in the rotor is turned by less than 9 °° in the meridian plane 
(see Fig. A(3c)). However, no data for p are available if this turning 
is 180 ° as in the present wheel, and for this reason the slip factor is 
taken as that of a standard centrifugal compressor wheel with e = 90 °. 

In accordance with Eq. A II ( 6 ) the flow coefficient cp, depends on p. 
For a chosen value p = 0.86 there is 

= (O. 86 )(cot 65 °) = 0.401 



For an assumed rotor efficiency = 0.86, 7 ^ = 17, and H^o/Rg = 0.55* by 
extrapolating Fig. A(3a) to higher values of cp^, there is p = 0.847. 
Hence, cp, = 0.395* and the slip factor for this velocity ratio is about 
O. 85 , giving a final value = 0.3964. For p = O. 85 , from Eq. A Il(5), 



W 2 

tT ^ = 0.6312 

VJ-, 

lo 



(3) 



Appendix A discusses the different possibilities to express the frictional 
losses in radial wheels. From Eq. A (25) 



c = 0.7200 

and from Eq. A(27) the so-called wheel efficiency T| w for T| R = 0.86 is 

Tl w = 0.4999 



With this value and Eq. 3, the velocity ratio 7 of Eq. A (2l) is 



W ? 

f - j r~‘ °-? 548 

2 is 

With these data the rotor dimensions at the discharge are obtained from 



Eq. A 11(11) 






(1.2739) 



(4) 



Hence, with Eq. 2 
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^2 ^2 1.2739 

R 2 ^ 11.7944 



0.1080 



The restriction factor at the inlet eye due to boundary layer growth 
is assumed to be O.98, and kg^ = 0.954, 

Thus 



iT = 0 • 110 9 (5) 

R 2 

In accordance with Fig. A(l6) this ratio is acceptable. 

The mean radius R^ is chosen as 5*5 in. , primarily to be able to 
drive the impeller by an available turbine and to mount the compressor 
to this item. Then the rotative speed N for = 800 ft/s is 

N = ^00 = (80011360) = l6j666 rpm 
tt( 5 . 5 ) 

30 12 



The basic compressor dimensions in inches are then: 



D 2 = 2 R 2 



D_ =21 
lo lo 



D. . =2 1. 
li li 



11.0 

6.05 

1.98 

0.6l 



The radial overlaps 6 of Fig. 1 at the impeller discharge are chosen to 
be 0.025 in., so that the blade height b^' at the wheel discharge equals 

b 2 ' = O.56 in. 

Drawing 2222 of Appendix B shows the design of the Hybrid compressor. 
Details of the impeller and the inducer are given in Dwgs . 2203 and 2304. 
The meridional channel of the wheel was laid out by assuming a blade tip 
contour with the smallest possible curvature changes between R^ q and the 
discharge. Drawings 2203 and 2204 show that such a contour could be 
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obtained with a two adjoining circular arcs with radii of 1.666 and 1.02 
inches. By assuming that the meridional flow areas, without the restric- 
tion due to the blades, change linearly with the distance X along the 
tip contour from at the inlet to ' at the impeller discharge the 
meridional channel could be laid out with the method explained in Fig. 2. 
As shown in Dwgs . 2203 and 2204 the hub contour thus obtained could be 
approximated by four adjoining circular arcs whose radii decrease from 
the rotor inlet to the rotor discharge. The clearance between the blade 
tips and the adjoining casing of Dwg. 2213 is everywhere equal to 0.025 
inch. 

2. Inducer 

Although some designers propose that the leading edges of the inducer 
blading have to be so arranged that particular incidence angles occur, it 
is shown in Ref. 2 that the correct incidence angles are obtained if the 
blockage because of the blade thickness is taken into account in the 
design. Figure 3 shows the conditions at the inducer entrance where V^, 
W_, are the flow properties just ahead of the inducer leading edges. 

In particular, the relative flow angle at the radius R^ q is equal to 
55° • At an arbitrary radius R^ the blade spacing is 



Because of the blade thickness t or its projected thickness t in peri- 
pheral direction, the axial component V 1 ' of the flow after entering the 
inducer is 




2 TT R. 



( 6 ) 




(7) 



u 



giving a flow angle 0^' of 



6 






= tan 



1 u) R . 



V 






(8) 



If the blade has this angle g^' at the leading edge there occurs an effec- 
tive positive incidence angle 

1 =■ h - h.' = h - % (9) 

Although in actuality the projected blade width is then t^' = t/cos g^' 
it is sufficiently accurate to take t^ = t/cos g^ because of the small 
values of i, and because the blade leading edge will actually be profiled 
as shown by the dotted contours of Fig. 3. For the same reason the thick- 
ness t is taken as 0.10 in., although the blade thickness t farther down- 
stream is actually 0.125 in.. With these assumptions the blade angle 
gL = g^' at the leading edge can be determined from Eqs. 6 to 8, or 

tZ R 



tan 



oo R_ t oo R., 

■i (1 - -?) - V- (1 - 



V. 



COS 3 1 2 TTR^ 



) 



Assuming to be constant along R^, and from the condition that at R^ q 

oo R 

tan 3 1q = 1 - = tan 55° 

there is with R_ = 3-025 in. 
lo 



tan ^ = g-i- tan 55° = tan 55° 



and 



tan gp = tan gJt — 1 
cos g 1 2 ttR 1 

Table I gives the angles g^, and the incidence angles i along the 
radius R^ at the inducer inlet. 

To be able to machine the impeller with simple means, it is assumed 
that its blade sufaces consist of generatrices that are everywhere parallel 
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to the radius at the inducer discharge where the blades must adjoin the 
radial impeller blades. As shown in Fig. 4 5 the camber line of the 
inducer blades can then be specified by a single function y = f(x). At 
a point P the blade surface has the angle y, where tan y = dy/dx. This 
angle , however , is the angle in planes A -A which are perpendicular to the 
radius at station where the inducer blades line up with the impeller 
blades. To obtain the actual blade angle Sg at a radius R^ the blade 
surface must be intersected with a cylinder of radius R^. For a point 
P(x>y) of the blade surface the blade angle 0^ at the radius R is from 

Fig. 4 



tan 



dy 



tan y 



cos 0 dx cos 9 



However 



cos 9 = 






IT 



K, 



Thus 



tan 



tan y 

B h - 



dy/dx 

JT- (y/Rp 2 



( 10 ) 



A function y = f(x) for the blade surface will now be established 
such that the blade angles Bg of Table I can be obtained along the radius 
by varying the axial length x of the inducer blades, as indicated in 
Fig. 3. In particular, the camber line of the inducer blades is chosen 
to be an ellipse. In accordance with Fig. 5 



y = b 






2 2 
a - x 



( 11 ) 



where 2a and 2b are the lengths of the major and minor axes of the ellipse, 
respectively. From Eq. 11 
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X 



( 12 ) 



^ = tan Y = - 

dx 1 a 



~rr 

V a - x 



2 1 



Referring to Fig. 3, the values of "a" and "b" will be determined by- 
specifying the conditions at stations C and D at the radii = 3.025 
in. and R^^ = 0.99 in. where, in accordance with Table I, the flow angles 
gg must be 50.32°, and 18.08°, respectively. As shown in Fig. 3, it is 
assumed also that the axial lengths x of the inducer blades are 1.8 in. 
at C and 1.255 in. at D. These data are sufficient to determine the 
values of "a" and "b". However if Eqs . 11 and 12 are introduced into 
Eq. 10, one obtains fourth-order equations which cannot be solved with 
simple means. For this reason a method of successive approximations 
is used by assuming values for the angles y T and of Fig. 5 at R^ q and 
R^ and calculating the corresponding angles gg by means of Eq. 10. Let, 
in accordance with Fig. 5, 

x T =1.8 = i 

Xjj = 1.255 = Qi = (0.697) i 
Then, from Eq. 12, with c = b/a, 



tan y T = 



tan v = 



(a /i) - 1 



.2 2 
C c 



“ (a /lf-f 



giving 



2 C^(tan 2 v m - tan 2 v„) 

I a\ b v *T T H 

w ' ~ 2 2 2 

C tan v T - tan t H 



(13) 






( 1 - C 2 ) tan 2 y T tan^ 



y H 



„2 . 2 

C tan V T 



tan 2 v H 



u*o 
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Successive approximations of y T and y H established the following final 
values : 

y T = 49.964°; tan y T = 1.1902 
Y h = 17.817°; tan y H = 0.3214 
For these angles, from Eqs. 13 and 14, 

7 = 1.04442; a = 1.8799 

JC 

J = 0.37460; b = 0.6743 

and, by Eq. 11 



^ = 0.37460 - 0.35867 J 1.09081 - (x/4) 2 ' (15) 

For x = l at R^ q = 3*025, 

y = 0.47973 

and for x = (0.697)/at R = 0.99j 

y = 0.1721 

With these values, from Eq. 10, at R^ q = 3-025, 



tan 



1.1902 

/l - ( 0.47973/3- 02 5 f 



1.20545 



and, at R = 0.99j 



tan p Bi = 0<32lU - , = 0.32637 

Jl - ( 0.1721/0. 99 ) 2 

Table I shows that these angles coincide very closely with the required 
blade angles at these radii. 

All blade angles fig of Table I at the different radii R^ can be 
obtained with the blade surface that is given by Eq. 15, for a particular 
change of the distance x of Figs. 3 and 5 with the radius R^. This 
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relationship between x and R^ must be determined with an iterative 
method. For chosen values of x there can be determined y and dy/dx 
from Eqs. 15 and 12. These values introduced into Eq. 10, for the radius 
R^ under consideration 
determined differs from the required blade angle of Table I, the distance 
x must be varied until the correct angle is obtained. Drawing 2204 gives 
the values of z = 1.8 - x vs . R^ that were obtained in this manner for 
the blade surface of Eq. 15 whose coordinates are listed also in 
Dwg . 2204 . 

With this design the inducer can be fabricated from a forging by 
first machining the tip controur of the blades and the surface of revolu- 
tion of the leading edges that is given by z vs. R. Then the blade 
passages can be produced by an end -mill that is moved in x and y directions, 
and along the radius, to produce the hub contour, However, with the 
described design the axis of the end -mill need not be tilted to produce 
the correct inlet blade angles. 

3. Diffusor 

The diffusor is arranged in a cylindrical annulus downstream of the 
impeller discharge, having the mean radius R^ = 5-5 in., and a radial 
height b^ = 0.6l inches. Ahead of the diffusor at station (2) the flow 
properties are obtained with the relations of Table A I. From Eq. A l(7), 
and the data established earlier, 

T 

= 1.08402 

x o 

From Eq. A I (8) 

T ' 

= 1.07013 

0 



, give the value of Eg* If the blade angle thus 
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From Eq. A I (6) 



From Eq. A I (20) 



and from Eq. A I ( 19 ) 



= 1.17411 

x 0 



TT = 1.2677 
0 



= 1.6764 
0 

The Mach number of the velocity V 2 ahead of the diffusor is from 
Eq. A 11(6) 

= 0.6446 

for the chosen flow angle 

a, = 65° 

A preliminary analysis showed that a minimum of two diffusor blade 
rows is necessary to produce an acceptable diffusor efficiency. These 
rows are designed in the same manner as bladings of axial-flow compressors 
by using the method of Ref. 3* In the following, station ( 3 ) is after 
the first diffusor blade row and station (4) is located after the second 
diffusor row, as shown in Fig. 1. 

Neglecting the change of the through-flow velocity because of the 
density increase from station (2) to station (3) the diffusion factor 
D of the first diffusor blade row is 



D 



cos ou (tan ql - tan a.) 

= 1 + 5 - 

cos 2 a 



cos 



°2 



(16) 



where a ^ is the average flow angle at station (3) and a the cascade 
solidity defined by 
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blade chord _ £ 
blade spacing - s 

For different deflections A a = and c it is possible therefore to 

evaluate D from Eq. l6. The losses increase with increasing values of D. 
For the first diffusor blade row these losses can be expressed by the 
total pressure loss coefficient defined by 



Y 



t 



AP _ P t2 - P t3 
P t2 - P 2 P t2 " P 2 



(17) 



where AP is the loss in total pressure in the blade row, and P^, P 2 j 
the total and static pressure, respectively, ahead of the cascade. Refer- 
ence 4 shows that Y of a stator blade row with a discharge angle oi^, 
diffusion factor D, and solidity <j can be calculated from 



Y = 



2 q 
cos 



°3 



0.0004 + 0.0639(D + O.l) 2 ' 91 + 0.057D 2 * 02 (1 - x) 3 ' 77 ] 



(18) 



The factor \ is defined by 



R -*H 



(19) 



where R is the radius at which Y is to be determined, and R_, R u are the 
tip and hub radii, respectively, of the blading. Since 1 - X = 1 at R^, 
and 1 - X = 0 at R^, the loss coefficient Y of Eq. l8 is higher at the 
hub than at the tip, if the other quantities in Eq. 18 do not vary radically 
along R. Equation 18 may give a realistic description of the loss variation 
in radial direction if the blade heights are large with respect to their 
chords, as for instance in axial compressors for jet engines. In such 
compressors the stator blades are attached at the tip and they usually have 
a clearance at the hub to prevent rubbing with the rotor assembly. The 
MSA calculating procedure of Refs. 3 or 4 does not separate the overall 
blade loss into profile losses, secondary or end losses, and tip clearance 
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losses. End and tip clearance losses depend on the ratio of blade height 

and chord and on the ratio of blade gap to blade height, respectively, 

and it is evident from Eq. 18 that these quantities do not influence the 

overall loss coefficient of the MSA design procedure. In the cascades 

to be determined the blade gap is zero but the ratio of blade height and 

chord, or bg/c, is quite small compared with the values in jet engine 

compressors, and the end losses will be very much larger than those in 

stator cascades of such machines. It is assumed therefore that Eq. 18 

for the mid-span, or for X = 0.5, establishes the profile loss coefficients 

Y only, or 
P 



Y = 2 a 

p cos 0^ 



0 . 00 k + o.o 639 (d + 0.1) 2 ' 91 + o.oo42 d : 



, 2.02 



( 20 ) 



To this profile loss coefficient will be added the so-called secondary 
loss coefficient Yg to give the total loss coefficient Y^, such that 



Y 



t 




+ Y„ 




( 21 ) 



Eq. 13(76) of Ref. 5, p. 380* gives the so-called induced drag coe- 
fficient of axial-flow cascades where the first part is the drag coe- 
fficient C—, due to secondary flow effects, whereas the second part is 
Db 

due to tip clearance losses . This part is taken to be zero because no 
radial blade gaps 6 exist in the present diffusor configuration. Hence 

C DS = °- 0k C L ° f < 22 > 

where the blade height h equals b^. The lift coefficient is equal to 
(see Eq. 13 (l4), Ref. 5, p. 337) 



2 

C T = — (tan ou - tan a 0 ) cos a 
La d 0 



(23) 



The angle a is obtained from 
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tan ou + tan ou 

tan = § (24 ) 

Equation 13(ll) of Ref. 5 relates the drag coefficient to the total 
pressure loss coefficient £ which, in the present application, would be 
defined by 



C = 



4 P 




It must be noted that £ relates the drop in total pressure to the discharge 
conditions of a cascade. Then, since 



Y s B P 



AP 



t2 ' P 2 



there is with Eqs . 22 , 23 , and 24 

Yg' = 0.l6 ^ cos 2 c &■> [ tan - tan o^] >/l + 4 [tan a, + tan o^] 2 ' ( 25 ) 

. ~ , 

l/cos a = + tan a m The secondary loss coefficient of Eq. 25 



since 



is denoted by Y_ * to distinguish it from a secondary loss coefficient Y " 
b £ 

that has been proposed in Ref. 6. In Ref. 6 it is stated that the ratio 

Y "/(Y + Y ") = Y "/Y," is a unique function of AQ , (c/h), where A a is 

s p S S "C 

the flow deflection in radians, and c/h the ratio of blade chord and blade 
height. This relationship is shown in Fig. 6. For known values of Y^ and 

Y "/Y." there are 
s ' t 



(Y "/Y ") 

v IT _ y ^ u 

s P 1 - (Y "/Y 11 ) 



(26) 



and 



v .1 P 

t - 1 - (Y//Y t ") 



(27) 



The losses in the diffusor blade rows will be determined with Eq. 25 as 
well as with Eq. 26 for purposes of comparison. These two methods are 
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based on different models and concepts. It can be noted that for given 
flow angles the secondary losses of Eq. 25 are primarily a function of 
the ratio of blade spacing and blade height, whereas with the method of 
Ref. 6, or Eq. 26, they are depending mainly on the ratio of blade chord 
and blade height. Evidently, if Eq. 25 is used, one obtains different 
cascade arrangements for a specified total loss than if Eq. 26 is applied, 
and this situation is indicative of the present state of ignorance about 
secondary flows in general. A survey of the available analytical methods 
for the determination of secondary flow effects is given in Ref. 7. 

As shown in Dwgs. 2202 and 2210-1, the diffusor blades will have a 
circular platform with a threaded cylindrical extension to attach them 
to the blade holders, which are shown as item 1 in Dwg. 2209. This design 
has been adopted to be able to change the blade stagger angle in the 
diffusor cascades if, because of the irregular flow conditions at the 
impeller discharge, the average flow angle differs from the design 
value of 65°. This diffusor blade attachment requires however, that each 
blade row does not have more than 38 blades. For Z = 38, the blade 
spacing s at the mean radius R^ = 5*5 in. becomes 

2 rr R 

s = — - — - = 0.9094 in. 

Zi 

and 

1 = — = 0 • 222it = 1.4908 
h b 2 0.61 

For particular values of a = c/s, 

£ = a 1 = ct(1.4908 ) 

n n 

For A a in degrees, the variable Aa(c/h) of Fig. 6 is, 



(28) 



(29) 



16 



o 



(30) 



Aa(c/h) = a oP a(l.^908) = a A a 

Table II shows the diffusion factor D, the profile losses Y , and the 
overall loss coefficients Y^' and Y^" as functions of solidity for flow 
deflections A a = 10, 12 and l4° in the first blade row of the diffusor. 

The quantities Y^' and Y^" are also plotted in Fig. 7. 

Although it is desirable to have the highest possible flow deflection 
A a in the first diffusor blade row, the general ' level of the losses as 
evidenced by Fig. 7 3 and the non-uniform and non-steady flow conditions 
at the impeller exit, make it advisable to limit A a to 12°. For this 
deflection the data of Fig. 7 show that the solidity should be about unity 
to obtain loss coefficients Y^ that are smaller than about 0.12. Experience 
with diffusors of centrifugal compressors seems to indicate however that, 
in spite of the results of Fig. 7 , more favorable conditions may be obtained 
for higher solidities, primarily because of the afore -mentioned irregular 
flow properties at the rotor discharge. At higher solidities the diffusor 
blades form channels with better guidance of the flow than if the solidity 
were low. These channels will act as flow straighteners which can equalize 
the flow irregularities with smaller flow separations than those which 
might occur in low-solidity cascades where the diffusion factors and the 
blade loadings are higher. It is likely then that the performance of 
the second diffusor blade row, and possibly that of the whole diffusor, is 
better than if the solidity of the first row is small. In general, the 
writer is of the opinion that the secondary losses do not increase as 
radically with the flow deflection A a and the solidity a as obtained by 
Ref. 6. 

For these reasons it was decided to build two sets of diffusors. In 
one, denoted by A, the solidity of the first row was chosen as 1.6. In 
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another, which is denoted by B, the first row has a solidity of O. 95 . The 
deflections A a have been chosen as 12° for both first rows of A and B. 
However, since the effects due to the non-steady impeller exit flows 
cannot be evaluated, the performance of these cascades will be determined 
with the loss coefficients of Fig. 7 • The design criteria of the two 
first diffusor blade rows are listed in Table III, which also gives the 
values of the chosen maximum blade thickness . 

The flow angle a ^ ahead of the second diffusor blade row is then 
65 - 12 = 53° • Theoretically, the highest pressure recovery would be 
obtained in the diffusor if this angle could be reduced to zero in the 
second row. It can be seen, however, that the blade loadings and the losses 
would become excessive for such deflections. Preliminary investigations 
showed that the deflections A a = 0 ^ - 0 ^ cannot exceed about 26° to 38 ° 
at reasonable loss coefficients, so that after the diffusor the flow angles 
will be between 27° and 15°. Even though the kinetic energy of these 
whirl components of the discharge flow at station (4) cannot be converted 
into pressure rise, the overall efficiency of the compressor will be better 
than the one obtainable if it were tried to reduce the flow angle 0 ^ to 
zero. 

Because it is intended to orient the blades of the second row in such 
a manner that they form a so-called tandem foil arrangement with the blades 
of the first row, the second diffusor row must also have 38 blades. Thus, 
its geometrical properties are likewise given by Eqs. 28, 29, and 30- The 
profile losses of the second row of diffusor blades can be determined 
by Eq. 20 if is replaced by 0 ^. The diffusion factors of these rows 

are obtained by Eq. l 6 , replacing by 0 ^ and a ^ by cv^. The overall loss 

coefficient is, similarly to Eq. 17, 
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( 31 ) 




The secondary flow loss coefficients Y will be determined by Eq. 26, 

s 

replacing by o^, and by - A a, and also with the method 

proposed in Ref. 6 by using the curve of Fig. 6. The results for different 
solidities and deflections A a are listed in Table IV and plotted in Fig. 8. 

Figure 8 shows that the overall losses increase very radically with 
solidity if the secondary flow losses are calculated in accordance with 
Ref. 6. If Eq. 25 is used however to determine Y , the overall losses Y, ’ 

S "C 

are almost constant for different values of cr, in fact they have a minimum 
for solidities of about unity. The data of Fig. 8 for Y that is, using 

- u 

Y " of Ref. 6, suggest that with solidities much smaller than unity a 
s 

diffusor with very low losses could be designed. However, not only would 
the blade chords be excessively small in such cascades , but the writer 
also believes that because of the high blade loading the flow would be 
too severely separated to make possible the small losses that Ref. 6 pre- 
dicts. If in contrast to this opinion, the method of Ref. 6 should prove 

to give the correct values of Y ", a cascade with solidity of 1.6, which 

s 

was chosen for blade row Al, would have excessive losses. It was decided, 
therefore, to design two second row diffusor cascades also. Both rows 
will have a deflection A & = “ % = 32 ° to keep the losses within 

acceptable limits, independent of whether the one or the other prediction 
method gives the correct loss coefficients. The so-called row A2 will be 
designed for a = 1.4, and for row B2 the solidity is chosen as unity. If 
Eq. 25 predicts the secondary losses properly, both cascades should have 
about equal losses. On the other hand, if Ref. 6 gives realistic secondary 
losses, row A2 would be greatly inferior to row B2. Although it will be 
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possible to arrange row A2 or B2 after either of the two first rows A1 
or Bl, for calculating purposes the so-called diffusor A is supposed to 
consist of rows A1 and A2 and diffusor B of rows Bl and B2. The design 
parameters of the second diffusor blade rows are listed in Table V. 

The blade profiles and their orientation in the cascade are determined 
with the method of Ref. 3- The necessary calculating steps are described 
in paragraph 8 of Ref. 8. The basic thickness distribution of the profiles 
is shown in Fig. 13 of Ref. 8. Except for slight modifications near the 
leading edge this distribution corresponds to that of a British C-4 pro- 
file. The camber line of the profiles is a circular arc with camber 
angle cp. For known values of cp and thickness ratios t/c the profile 
coordinates are obtained with the relations of Fig. 14 of Ref. 8. 

Figure 9 shows the attitude of a blade in a diffusor cascade and 
explains the symbols that are used for the profile calculations in Table 
VI. By using the method of Ref. 3 one often obtains excessively large 
incidence angles. As seen from Table VI this situation occurs for row 
A1 where i = + 4.312°, for row Bl where i = - 2.95, and for row B2 where 
i = - 7.002' . For these rows, incidence angles i* of +2°, -1°, and -2°, 
respectively, are chosen, and the corresponding camber angles cp* and the 
stagger angles y* are calculated with the values of d6/di of Fig. 177 
of Ref. 3* Table VI also lists the drawings of Appendix B where the 
different blade profiles and their coordinates are shown. From these 
drawings the angles Ay of Fig. 9 were determined graphically and were 
used to establish the stagger angles y' that are necessary for the installa- 
tion of the blades in the blade holders of Dwg. 2209. 
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II. DESIGN POINT PERFORMANCE 



1. Rotor 

For a chosen total temperature Tq = 520°R at the compressor inlet, 
and for U^/a q = 0.7156, = 0.4842 (see Eq. l), the static temperature 

T^, and the static pressure at the inlet eye of the impeller are 
obtained from Eqs. A I (5) and A I (l8 ) ? or 

T 1 

ijr = 0.9848 
0 

P n 

TT = 0.9478 
0 

The absolute value of the total pressure P^ at the impeller inlet is as 

yet 'unknown because in the test rig the compressor operates as an exhauster. 

Drawing 2223 shows that the air entering the compressor is expanded in a 
throttle valve to a pressure below atmospheric pressure. From a plenum 
chamber (Dwg. 2164-1 ), where the kinetic energy of the flow through the 
throttle valve openings is destroyed by a system of adjustable screens , 
the air enters the suction pipe through a profiled orifice. This suction 

pipe with an inner diameter of 18 inches has a length of about 18 feet, 

so that non-uniformities of the flow which have not been removed in the 
plenum chamber will be equalized ahead of the flow measuring nozzle of 
Dwg. 22l6. This nozzle has a diameter of 5-38 inches, hence its flow 
area is only about 9 percent of the cross sectional area of the suction 
pipe. The large acceleration of the flow in the nozzle will produce a 
uniform flow at its discharge. Drawing 2222 shows the location of the 
flow nozzle with respect to the compressor inlet. Drawing 22§2 shows that 
provisions have been made in the nozzle casing (Dwg. 2l4o) to carry out 
flow surveys ahead of the 7 in. diameter duct ahead of the impeller by 
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means of Pitot probes. 



Drawing 2222 shows that a honeycomb flow straightener is arranged 
after the diffusor blade rows to remove the whirl components of the flow 
at this station, and that the air is then discharged into the atmosphere. 
This arrangement has been adopted for simplicity of the design, since it 
does not require a scroll or a complicated collector after the compressor 
from where the flow is discharged into the atmosphere through a throttling 
valve. Moreover, with this arrangement the driving power of the compressor 
is lower than with a test set-up where the inlet flow of the compressor 
is at atmospheric pressure. However, the total inlet pressure P^, which 
is equal to the total pressure in the 7 in. diameter inlet duct, depends 
on the pressure ratio that is produced by the compressor. If, as shown 
in Fig. 1, the conditions at the exit of the discharge duct, arranged 
after the annular flow straightener, are denoted by the subscript ( 6 ), 
the static pressure pg must be equal to the ambient atmospheric pressure 



P , , or 
atm 



P = P 
0 atm 



(p 6 /p o ) 



( 32 ) 



Hence it is necessary to determine the pressure ratio P^/Pq to obtain P^. 
Evidently, P^/Pq depends on the pressure drops in rotor, diffusor and 
flow straightener which will be evaluated in Section IE. 3. 

The weight flow rate w through the compressor can be determined from 
Eq. 2. For the assumed blockage factor kg^ = O. 98 , and since Rg = 5.5 
inches , 



9 JL 

p o 



IE„ T„ TT Rg kg 1 



G '“0 _ ' “2 J1 B1 _ „ o_,_ ( . 2x 

1 t jj cy | )| 7*8963 (m. ) 



( 33 ) 



With R g = 53.35 ( ft -lb ) / ( lbm , °R ) , g = 32*174 ft/s 2 , T Q = 520° R, 
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w = 0.2689 P Q (lbm/s ) (34) 

In Eqs. 33 and 34 , Pq is the absolute total inlet pressure in psia. 

In the rotor the total temperature of the fluid is increased from T Q 
to T^2 • an adiabatic process the energy increase due to this tempera- 

ture rise is equal to the energy per unit flow rate that the rotor must 
transmit to the fluid. Hence, the power HP required to drive the rotor 
is, exclusive of bearing and other mechanical losses, 

" c p (T t2 - T o> 5?0 (HP) (34a) 

With 



c = R, 



p “G Y - 1 



(ft-lb)/ (lbm, °R) 



and Eq. 33 



HP 



- 7.8963 Y p r ~ R - T v ^t2 , 
550 v - 1 F 0 * S R G T 0 



550 Y - - ^ — 0 

For y = 1.4, and since T^/Tq = 1.17411, 



HP = 8,7489 



P Q V S r g t 0 (10 3 ) 

For T Q = 520° R, y = 1.4, R p = 53*35> there is also, 



G 



HP = 8.2656 P r 



(HP) 



(35) 



(36) 



In Eqs. 35 and 36 the pressure P Q must be in psia. 

The following data at the rotor discharge were obtained earlier: 



= 1.08402 



= 1.2677 



t2 

■0 



= 1.6764 
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“V2 * 4 ' °' 6W 

With T_ = 520°R, or a Q = 1117 . 85 ft/s, the velocity of sound a 2 at the 
rotor discharge is 



a 2 = a Q J T 2 /T 0 ' = II 63.86 ft/s 

and 

V 2 = 750.23 ft/s 

This average velocity at the rotor exit has an average flow angle = 6 ?°. 
2 . D iff us or 



The total pressure P._ 
from Eq. 21, or 



after the first diffusor blade row is obtained 




(1 - 



Y ) — 
^ P 0 



Y 1* 

* p o 



(37) 



where the pressure loss coefficients Y^ of the rows A1 and B1 are listed 
in Table III. At station ( 3)5 after the first diffusor row, the velocity 
has an average flow angle = 53° • Although the annulur flow areas 
are equal at the station (2) and ( 3 ), the boundary layer growth in the 
blade row will make the effective flow area smaller at ( 3 ) than at (2). 

At station (2) a restriction factor k^ 2 = 0.954 has been assumed earlier 
to establish the diffusor annulus, or 

A 2 - ^B2 A geom 

where A is the cross-sectional area of the annulus. At station (3) 
geom 

it is assumed that the effective flow area A^ is 

A 3 = *B3 A geom 

and k _ 0 = O .92 by assumption. 
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The flow conditions at station (3) must be determined analytically 



for the flow area A^, the flow angle o^, and the total pressure P^, by 
using the energy equation, the equation of motion, and the equation of 
continuity. The energy equation for an assumed adiabatic process of a 
steady flow gives 




In fact, the same total temperature will occur at stations (4), ( 5 ) and 

(6), if the flow processes are adiabatic. 

Since the conditions after the second diffusor row at station (4) 

must be established with the same approach as those at station (3), a 

calculating procedure will be developed with the symbols of Fig. 10. 

Figure 10a shows a diffusor cascade in an annulus. The effective flow 

areas are A = k_ A and A, = k_, A at entrance and discharge, 
e be geom d Bd geom 

respectively. The average axial components V eg and V^ a that are uniform 
over A g and A^, respectively are 



V 


= V 


cos 


a 


ea 


e 




e 


V da 


>* 

1! 


cos 


a d 



A section through the blading of the diffusor row is shown in Fig. 10b, 
and Fig. 10c depicts the thermodynamic process between entrance and discharge 
in an entropy diagram. From the equation of motion 



T,. = T., = T + 



te td e 2 g c 



= T, + 



d 2 g c 



P P 

From the equation of continuity, with p = p/(g E„ T) 



(38) 



A d V d C ° S «dR G T d 



P d P e 

a = A V cos a 



e e 



e R„ T 
G e 



With A^/A^ - Eq. 39 can be rewritten as 



(39) 
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V, 



/« V E g T a ’ 



V e ,/< V T q ) ' 

‘‘Bd C ° S “d J (T e /T 0 )' VsvRj T e * P d /P 0 ) 



(UO) 



With the Mach numbers 



> 

11 

2* 


(41) 


Jg Y R g T d ' 




V 




M = 

e YivR G T e - 


(42) 



and 



cos oi 
e e 



^ 

k B4 C0S 



V P C 



“d J (T /T.) 



M 



(43) 



which is a constant for known entrance conditions and area ratios, there 
is from Eq. 40, 



Md = K 



J < V T 0> ' 
W 



From Fig. 10c at the entropy s^, from Eq. 38 , 



(44) 



and 





(1 + 




1 ) 



From these relations 




T 



1 + 



t Jl 

Y - 

2 



0 

1 




(45) 
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and 




(46) 



where, similar to Eq. 37 > with the pressure loss coefficient Y of the 

"C 

blade row, 



P 



td 




(47) 



Because of their exponential expressions it is not possible to substitute 
Eqs. 45 and 46 Into Eq. 44 to obtain a closed solution for M^. It is 
necessary to choose values of M^, say M^*, to calculate T^/T^ and P^/Pq 
from Eqs. 45 and 46. These quantities are then used in Eq. 44 to deter- 
mine the Mach number M, . If differs from M *, other values of M * 

d d d d 

must be chosen until agreement is reached. With the final value of 

M, = M * it is possible to calculate the ratios T^/T^ and p,/P„ with 
d d d' 0 d' 0 

Eqs. 45 and 46, and the velocity is obtained from Eq. 4l, or 



This iterative calculating process has been carried out for the rows A1 

and Bl, for both overall loss coefficients Y^ 1 and Y^" that are listed 

in Table III. In Table VII where the final data are presented, the 

columns designated by Al* and Al" give the performances of blade row A1 

for the overall loss coefficients Y^' and Y^", respectively, of Table III. 

A similar designation is used for the performance of blade row Bl, for 

the losses Y ' and Y, ". With the flow conditions thus determined at (3), 

"C *c 

the flow properties at station (4), after the second diffusor blade row, 

can be calculated with the same method. These calculations are also 

carried out with both loss coefficients Y ' and Y, " of Table V. The 

t t 




(48) 
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quantity K of Eq. 43 is determined for k„ 

Be 

The flow angles and are 53° and 21°, 



= kj, 3 - 0.92 and = 

respectively, in accordance 



0 . 88 . 



with Table VI. 



The calculations of the flow conditions after the diffusor blade rows 
are quite lengthy and have to be carried out with great precision. On the 
other hand they do not warrant the use of a high-speed digital computer 
for which a rather elaborate program would have to be set up to perform 
the necessary iterations automatically. For the debugging of such a 
program several tries would probably be necessary until results could be 
obtained. However, for evaluations of the type discussed above, in fact 
for most calculations of this report, modern electronic and programmable 
display calculators are admirably suited. The Monroe Model 1655 with Card 
Reader CR-1 used at the Turbo- Propuls ion Laboratory has proved to be an 
extremely valuable tool, and the writer believes that such calculators 
will completely revolutionize analytical engineering approaches. Libraries 
of programs for a variety of engineering calculations are already available, 
and programs for particular tasks can be established with great ease. 

Since most calculators, specifically the Monroe Model 1655, have provisions 
to verify and debug programs with a step by step procedure, program errors 
can be detected and corrected without delay by the programmer himself, 
without having to go to the trouble of submitting several modified 
programs as is necessary in most central computer facilities. At present, 
many high-speed computers are used for programs that do not really need 
the small access time and the large storage capacity of modern systems, 
and the use of programmable calculators will make available more time for 
those calculations that can be solved only by high-speed machines. Con- 
versely the engineer or student is relieved of many routine calculations, 
without having to worry about inaccuracies and errors, and he obtains 
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results when needed during the progress of his work. Thus more time is 
available for creative efforts where the formulation of a problem is more 
important than the manipulation of numbers that is necessary for its solution. 

To show the possibilities that exist in solving engineering problems 
with the Monroe 1655 calculator, some of the programs used for the calcu- 
lations for this report are enclosed in Appendix D. Each program consists 
of a description, an operating instruction, and a listing of the program 
steps with their codes. The serial numbers of these programs refer to 
the private program library of the writer. The data of Table VII for the 
conditions after the two diffusor blade rows have been obtained with 
Program Wo. 109. 

Appendix D also contains a more elaborate calculating program that 
requires the use of the Monroe calculator Model 1880 which is an advanced 
and more expensive machine with larger memory capacity. 

The pressure ratios P^/Pq of Table VII will be used to determine the 
so-called total-to-total compressor efficiency T^. For an isentropic 
compression from the total inlet pressure Pq and the total inlet tempera- 
ture Tq to the total discharge pressure P^, the rise in total temperature 
AT^ equals 



* T tis ■ T 0 [<W (Y ' l!/Y - l] 



Since this temperature rise is proportional to the theoretical energy in- 
put per unit mass flow rate, and because for the actual process the energy 
per unit mass flow rate necessary to drive the rotor is proportional to 
T^ - Tq, which equals T ^ - T^, the compressor efficiency is 

<VV <Y ' 1)/Y - ^ 



\4 = 



< w - 1 



(49) 
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The compressor efficiencies 7^ are given also in Table VII. If the 
secondary flow losses are calculated in accordance with Eq. 25, an effi- 
ciency of 81.5 $ is obtained for the compressor with diffusor A, whereas 
if the losses for this configuration are determined with the method of 
Ref. 6 the efficiency is 78 . 9 $ only. The calculated compressor efficiencies 
with diffusor B are 82.3 $ and 82.8$, respectively, with the two methods. 

It would appear, therefore, that higher efficiencies can be reached 
with diffusor B than with diffusor A. However, for the reason mentioned 
earlier, it may be that more favorable conditions can be obtained with 
the high-solidity diffusor A. Moreover, the performance calculations 
cannot take account of the radial velocity gradients or the non-steady 
flow at the rotor discharge, and it is because of these uncertainties 
that tests must be carried out to verify which diffusor configuration 
gives the best performance. 

3 . Flow Straightener and Discharge Duct 

The flow conditions in the flow straightener and the discharge duct 
of Fig. 11a can be determined with approximate methods since the Mach 
numbers of the flow are small. At station (4), at the diffusor discharge, 
the velocity has the angle = 21° with the axial direction. The 

peripheral component sin will be destroyed as the flow enters the 
axial flow straightener. This element consists of a honeycomb with hexa- 
gonal channels that have a spacing of 0.125 inch, a wall thickness of about 
0.002 to 0.003 inch, and an axial length of I .38 inch. As shown in Dwg. 
2212, twelve honeycomb segments are arranged in an annulus between spacers 
that are 0.25 inches thick. On assuming that the flow leaving the diffusor 
does not diverge in radial direction, the velocity Vg in the channels of 
the flow straightener is with the symbols of Fig. 11a, neglecting density 
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changes , 



\ 2 nK,, b 2 = V s [2 ttR, t> 2 - (12)(0.25) b,, ] (§fj§|) 

where the fraction 0.123/0.125 is the blockage due to the wall thickness 
of the honeycomb straightener . With V^ a = cos 21° = 0.933 V^, R^ = 5-5, 



V . y, cos 21° 1 2 > 1 ?? 10 L 123 ) . 1-039 
S 4 , (12) (0.25) 4 

2 tt (5.5) 



(50) 



As shown in Fig. 11a, it is supposed that after the honeycomb the flow 
will diverge radially to a height b,_ = O .65 in. at station ( 5 ). Hence 
the velocity is assumed to be 



v 5 - bf V * &§ (°-933) v 4 - 0.8764 v k 



(51) 



Because the discharge duct between stations (5) and (6) is not a properly 

designed diffusor it will be assumed that the flow will follow the dash- 

dotted boundaries of Fig. 11a. Neglecting the change in density from (5) 

to (6), the velocity is then 

0 

v 6 = v 5 I " - °- 824 v 5 - °- 722 \ 

Figure lib shows the flow process between stations (4) and (6) in an 
entropy diagram. If the flow is adiabatic the total temperature is 
everywhere constant, equal to T^ = T^. From the diffusor discharge at 
(4) to the entrance of the honeycomb the velocity is reduced to V 4a at 
the constant static pressure p^. This process is associated with an in- 
crease in entropy from s^ to s^', reducing the total pressure from to 

P t4’* The velocity is then increased from V^ a to at station (4') 

just inside the entrance to the flow straightener. If this acceleration 
is isentropic the pressure p^ is reduced to p^'. With the symbols of 
Fig. lib, 
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2 2 

* \ - V 4a 

T 4 =T 4 + — ~ 

2 g c 



Since = R Q y/(v - l), and V, lo = V,, cos o), , 



4a 



V 



where = 



V = \ + ^JL A sin °4 \ 

T o T o 2 ' a o 

1117.85 ft/s and =21°. Further, 



P t4* P 4 /^t2^ T 0 



y/(y - 1) 



p o F o v > J -4" /x o / 



/_t2___0\ 

vv7v 



The velocity increase from to Vg at the entropy s^' reduces the 
temperature from T^* to ' , where 



V„ 



T 4* = T t2 



2 g c 



or, with Eq. 50, 



Also 



and 



V T 

r o ^0 



Y - 1 A° 39 A 
T T 2 V a 0 / 



V 



t4 



r “( 4 ) 



, v/ (v - 1) 



V V A'/T, 



0 



r 4 /x o \ 

P 0 ' P t2^ T 0 ' 



y/(y - 1) 



The loss in total pressure A Pg = P^' - P^" in the honeycomb will 
determined with the frictional coefficient f of flows in pipes, or 



( 53 ) 



(54) 

static 



(55) 



(56) 

be 
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( 57 ) 



A P = (— ) f — V 2 

a s Sa^ 2 s 



The hydraulic diameter of the flow channels in the honeycomb 
straisiitener is 0.123 in., and L = I .38 in. The factor f depends on the 
Ticynolds number 



R = 
e 






Table VII shows that the velocities do not differ greatly, hence R g 

can be determined for an average value = 330 ft/s. Moreover the average 

ratio Tj^/Tq is about 1 . 157 ? and the pressure p^ is nearly equal to the 

atmospheric pressure. For T^ = (1.157)(520) = 602° R, and p^ = lU.7 psia, 

h 0 

there is v = 1.97(l0~ ) ft/s . Hence, with Eq. 50, 

H (]-039)(330)(0.l23/l2) _ 1 . 7 o (1 ^ ) 
e 1.97(lo' 4 ) 

From Ref. 9 for smooth surfaces, f = 0.0265, and by Eq. 57 

4 P S - (0 - 0265) i v s 2 ' °- 3 ? v s 2 

At the low Mach numbers of the flow, there is approximately 



1 V S 2 = P t4' - p 4* 



or 



V V p 4* 
T" = °-7 -Ir- + °- 3 T" 

0 0 0 



(58) 



From Fig. lib, and Eq. 56 the static pressure p^" is then obtained from 



V / p t4"/ p oV p 4\ 

p o 



( 59 ) 
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The reduction of the velocity from V G to V is supposed to occur at the 
static pressure p^", hence p,_ = p^", as indicated in Fig. lib . The static 
temperature is with Eq. 51, 



t2 



and the pressure ratio P /P 




^0.8764 v 4 2 

V ^ ) 



is 



(60) 



P + * P)," / T + o/T 



t5 _ ( L t2 /x o\ 

p o ' p o WW 



v/(v - l) 



(6i) 



o ^0 ' • l 5 / o 

The loss in the discharge duct between stations ( 5 ) and (6) is assumed to 
be 

P t6 - P t5 = <°- 4 >< P t5 - V 

Thus, with p = p^" 



0.6^ + 0.k P ± 

0 0 0 



With Eq. 52 



(62) 



T 6 T t2 Y - 1 v 4\ 

T 0 “T 0 " 2 V a. ) 



0.722 V, 



and 



n P T /T ^/(y “ -4 ) 

= _t6 / t 6 /t o \ 

p o~ p o 



(63) 



(64) 



For the data at the diffusor discharge (station (4)) of Table VII the 
conditions of state in the flow straightener and the discharge duct have 
been calculated with the above relations. The results are summarized in 
Table VIII. As discussed earlier, the static pressure pg at the exit of 
the duct will equal the atmospheric pressure. If this ambient pressure 
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P atm '*' S taken as psia, the total pressure P Q at the compressor inlet 

is obtained from Eq. 32. The last line of Table VIII gives these pressures 
for diffusors A and B and the two different losses which were taken into 
account for their performance prediction. 



4. Compressor Performance at Design Point 

The weight flow rate w and the driving power of the compressor can 
be determined from Eqs . 34 and 36 for the values of P^ from Table VIII. 
The pressures thus obtained are for a total, inlet temperature T^ = 520° R 
and a speed of rotation N = 16,666 rpm. Of interest is also the torque M 
that is necessary to turn the compressor, which equals 



= HP(550) = HP(550) 

0) TT N/30 

With Eq. 34, and N = 16,666 rpm 



(ft-lb) 



(65) 



M = 2. 6045 P Q (ft-lb) 



or 

M = 31.2554 P Q (in -lb) ( 65 ) 

where P^ must be in psia. In Table IX the compressor performance data 
for the design point at N = 16,666 rpm and T^ = 520° R are summarized. 

Also given are the pressure ratios P^/P^ and p^/P^, and the efficiencies 
with the two diffusors A and B, for the two methods used for the 
determination of their secondary flow losses. 

It is evident from Appendix A (Ref. l), and the calculations of this 
report, that t^/a^, y = c^/c^, and the flow angles are the prime variables 
that establish the performance of a particular compressor, in fact of all 
designs which are geometrically similar. Except for possible changes in 
the loss coefficients which will be discussed later, the pressure ratio 
P^/Pq of these machines will be equal if U^/a^, y, and the flow angles 



35 



are equal. This situation occurs because the density ratios of the fluid 
and the Mach numbers of the flow remain constant at particular stations 
in the compressors, which makes it possible to maintain the sarnie flow 
angles with respect to the bladings. From Eq. A Il(9) it is seen that 
the dimensionless mass flow rate (m)* defined by 



(m)* 





(66) 



must be constant at these conditions. From Eq. 35 and Eq. A 1(6) it is 
apparent that the dimensionless power (HP)*, defined by 



(HP). = ( HP > (55°) 

n R 2 P 0 R g V 



(67) 



is invariant also for similar machines that operate at the same values 
U^/a^ and y, and produce the same pressure ratio P ^/P . 

Since uo = the moment of Eq. 65 becomes with Eq. 67 , and 



a 0 g T q T q ,' 



M = 




(HP)* 

(U 2 / a 0 ) /7 



( 68 ) 



Hence, similar to Eqs . 66 and 67 the so-called dimensionless moment (M)* 
will be defined as 



(M)* ^ (69) 

"V P 0 

which is constant also for the conditions at which (m)* and (HP)* do not 
change . 

Obviously the dimensionless quantities of Eqs. 66, 67 , and 69 have 
equal magnitudes in any consistent system of units for geometrically 
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How- 



similar compressors that operate at the same values of Ug/a^ , y, and 



P t4 /P 0 , independent of their size and the magnitudes of P Q and T Q . 
ever, the losses in these machines do not depend only on the angles and the 
Mach numbers hut also on the Reynolds numbers of the flow. That the Mach 
numbers remain constant at fixed values U^/a^, y, and flow angles, indepen- 
dent of the magnitudes of T^ and P^, is shown by Eqs . A 11(4) and A Il(6). 
The flow angles, and therefore the blade incidence angles, remain equal 
at particular stations because the ratios T/T ^ and P/P^, which establish 
the density ratios p/ p^, are also not affected by Tq and Pq. The Reynolds 
number R g at a particular station, where the velocity V exists at p and T, 
can be expressed by 




since, for geometrically similar compressors, the dimensions of all flow 
channels are constant fractions of the mean compressor radius R^. If M^ 

is the Mach number of V, or V = M^ J y g R^ T, and, since p = p/(g R^ T), 
there is also 



of units used in this report, the term in the square bracket is dimensionless. 
This is necessary because the term in the round bracket and R g are dimension- 
less also. 

Whereas, for reasons discussed earlier, the expression in the round 
bracket is constant for all geometrically similar compressors operating at 
particular values of U^/ a^ , P^/Pq and y, independent of P^ and Tq, the 
term in the square bracket depends on the compressor size, the gas properties, 




( 70 ) 



2 

Since the dynamic viscosity y, has the dimension (lb-s)/ft in the system 
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and the magnitudes of P Q and T^. Experience has shown, however, that the 
influence of Reynolds number on losses in turbomachines is quite small, 
and that large increases of the losses occur only below some critical value 
of R g which is usually reached only at very low pressures or with small 
dimensions. For most gases the dynamic viscosities increase slightly 
with temperature but the absolute magnitudes of ^ do not vary greatly 
for different gases. Most diatomic gases such as N^, 0 Hg, CO, and 
also air, have values of y of about 1.4 at ambient temperature, and the 
ratio of the viscosities of the heaviest of these gases (Og) and the 
lightest one (H^) is only about 2. A change in Reynolds number by a 
factor of two has a very small influence on the losses, except near the 
above-mentioned critical value of R . 



This discussion shows that geometrically similar compressors that 
operate with gases that have the same specific heat ratio will have almost 
equal values of (m)*, (HP)*, and (M)* at particular values of U^/a^ and 
V P 0> if the Reynolds numbers are sufficiently high. Evidently the 
efficiencies of the compressors at these operating conditions will then 
also be nearly equal, irrespective of the molecular weight of the gas and 
the magnitudes of inlet pressure and inlet temperature. The quantilUS (m)*, 
(HP)* and (M)* are shown in Table IX for the ratio U^/a^ = 0.7156. 

For a particular compressor operating with a particular gas, the 
dimensionless quantities of Eqs . 66, 67, and 69 can be simplified by 
introducing the ratios 



T 



®C T. 



0 



REF 



and 



(71) 



6 C = 



0 



REF 



(72) 
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where T and P are usually taken as 518.4° R and 14.7 psia. The 

KHjI 1 KHjr 

subscript C was introduced to differ ientiate these ratios from similar 
ones that later on are used for the performance calculations of the drive 
turbine. 

With Eq. 71) the ratio Ug/a^ is also 

U 2 (tt/30) N R 2 

a o / ’ 

7 g y r q e c 

Hence in a compressor that handles a specific gas, equal values of U^/a^ 
are obtained for equal ratios N/7""9^» Thus, 

= N//0“ (rpm) (73) 

will be called the referred speed of the compressor. Similarly, from Eq. 
66, the so-called referred weight flow rate is defined by 

W REF = bl (lbm/s) (74) 

u 

From Eq. 67 the referred power is 

(HP)^ = — (HP) (75) 

6 C 7^ 

and from Eq. 69 the so-called referred torque becomes 

M 

Mppp = — (ft-lb) or (in-lb) ( 76 ) 

6 C 

On the basis of the previous discussion it can be stated that a compressor 
designed for and operated with a particular fluid will have the same 
pressure ratio at the same values of IL,-,- and w independent of the 
absolute magnitudes of P^ and Tq, provided that the possible variations 
of Reynolds number have a minor effect on the losses in the bladings of 
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the machine. Then the referred power required by the compressor will 
also be a unique function of N /J 8 C and w J 9q/ 6 c , indicating that its 
efficiency is only depending on these two variables also. The main dis- 
advantage of the so-called referred parameters is that they are not 
dimensionless, and that they cannot be used to compare the operating 
performance of different machines with each other, as it is possible with 
the earlier defined non-dimensional quantities . 

The referred operating parameters of Eqs . 73 to 76 for the design 
point are given in Table IX. 
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III. OFF -DESIGN POINT PERFORMANCE 



1. General 

Because the Hybrid compressor operates as an exhauster in the test 
rig it is of interest to know its performance at pressure ratios and 
speed other than those for the design point. If the pressure ratio is 
smaller than the design value, the pressure P^ will be higher and since 
the weight flow rate, and the driving power are proportional to P^, it is 
likely that the maximum power occurs at pressure ratios lower than the 
values at the design point. 

2. Rotor Performance 

Prime variables for the rotor performance are the peripheral speed 
ratio U 2 /a Q and the absolute discharge angle a^. If particular values of 
these quantities are chosen and if the slip factor p, and the velocity 
ratio y were known it would be possible to determine the dimensionless 
mass flow rate through the wheel by means of Eq. A Il(ll). These calcu- 
lations can be carried out by means of program 102 of Appendix D, which 
also establishes the ratios T^/T^ and p^/P^ for the assumed values of p, 
and y • Actually, program 102 determines the quantities: 



and 



w J(Rjg) T 



0 



= (m) tt R, 



‘0 



W V(R r /g) T n w V(Rp/g) T 



= 



0 



*1 = 



2 A 2 ^B2 p 0 2 tt R 2 b 2 ^2 P 0 

w V(R G /g) T Q 



w V( R G /g T 0 



A 1 *B1 p o 



" ( R 10 2 - R li 2 ) ^Bl p 0 



(77) 



(78) 



(79) 



for chosen blockage factors k^ and kg 2 . 
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The dimensionless flow function of Eq. 79 can be used to estab- 
lish the ratio of the static pressure and P^, since is identical 
with the total pressure P at station (l), with the relation 

/ o f 2/ y (y + l)/yi 

'r^lw - W ] < 81 

Since for known values of and y, it is not possible to calculate 
P 1 /P 0 in closed form, the calculating program 103 of Ref. D has been set 
up to determine p^/P^ with an iterative process. 

Program 104 of Appendix D can then be used to determine the flow 
properties at the inducer inlet with the relations of Appendix A which 
are listed in the program description. This program gives T^Tq, V^/ a Q , 
Wio/ a o> and the quantities and i'. The latter is the incidence angle 
at the outer inlet radius R^ q if blockage due to inducer blade thickness 
is taken into account. 

Appendix C analyzes available test data of losses in centrifugal 
compressor rotors, and shows that Fig. 24 can be used for the present 
off-design calculations. From Fig. 24 the rotor efficiencies are 
obtained as functions of i' and M^. The efficiency is defined by 
Eq. A (12). With Eqs. A Il(5), A (27), and A (21) it is possible to 
calculate the velocity ratio y = W^/W^^ which is used for the perform- 
ance calculations with the method of Appendix A. This conversion can be 
carried out with program 105 of Appendix D. The wheel efficiency T]^ of 
Eq. A (15) and the deceleration ratio W^/W^ are outputs of this program 
also. The value of y thus obtained should coincide with that initially 
assumed for use in program 102, and y must be iterated by successive 
approximations till agreement is reached. 

For the design point the slip factor p, = jj, = 0.85 has been deter- 
mined from Fig. A (3a) for the given rotor dimensions, the assumed rotor 
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